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New results after 0.69 kgyr of measurement with an enriched 76 Ge detector of the Heidelberg- 
Moscow experiment with an active mass of 2.758 kg are presented. An energy threshold of 9 keV 
and a background level of 0.042 counts/(kgdkeV) in the energy region between 15 keV and 40 keV 
was reached. The derived limits on the WIMP-nucleon cross section are the most stringent limits on 
spin-independent interactions obtained to date by using essentially raw data without background 
subtraction. PACS number(s): 95.35. +d, 14.80.Ly 



INTRODUCTION 

The nature of dark matter in the Universe remains 
a challenging question. Even if new measurements will 
confirm that we live in a low f2 ma tter universe (Omatter ~ 
0.3-0.4) |,| a considerable amount of nonbaryonic dark 
matter is needed. WIMPs (weakly interacting massive 
particles) are among the most discussed candidates j|, 
being well motivated from early universe physics Q and 
supersymmetry ||. 

WIMP detection experiments can decide whether 
WIMPs dominate the halo of our Galaxy. For this reason, 
considerable effort is made towards direct WIMP search 
experiments which look for energy depositions from elas- 
tic WIMP-nucleus scattering ||. Germanium experi- 
ments designed for the search for neutrinoless double beta 
decay were among the first to set such kind of limits [frfi . 

The Heidelberg-Moscow experiment gave the most 
stringent upper limits on spin-independent WIMP in- 
teractions j|] until recently. The present best limits on 
the WIMP-nucleon cross section come from the DAMA 
Nal Experiment jT3|. The Heidelberg-Moscow experi- 
ment operates five enriched 76 Ge detectors with an active 
mass of 10.96 kg in the Gran Sasso Underground Labo- 
ratory. It is optimized for the search for the neutrinoless 
double beta decay of 76 Ge in the energy region of 2038 
keV. For a detailed description of the experiment and 
latest results see [fl2f . 

In this paper we report on results from one of the 
enriched Ge detectors, which took data in a period of 
0.249 years in a special configuration developed for low 
energy measurements. A lower energy threshold and a 
smaller background counting rate has been achieved with 
the same detector as used in 1994 |9), mainly due to the 
lower cosmogenic activities in the Ge crystal and in the 
surrounding copper after four years without activation. 



EXPERIMENTAL SETUP 

The utilized detector is a coaxial, intrinsic p-type 
HPGe detector with an active mass of 2.758 kg. The 
enrichment in 76 Ge is 86%. The sensitivity to spin- 
dependent interactions becomes negligible, since 73 Ge, 
the only stable Ge isotope with nonzero spin, is deen- 
riched to 0.12% (7.8% for natural Ge). The detector has 
been in the Gran Sasso Underground Laboratory since 
September 1991; a detailed description of its background 
can be found in @. 

The data acquisition system allows an event-by-event 
sampling and pulse shape measurements. The energy 
output of the preamplifier is divided and amplified with 
two different shaping time constants, 2 /xs and 4 /zs. The 
fast 2 fis signal serves as a stop signal for the 250 MHz 
flash analogue to digital converter (FADC) which records 
the pulse shape of each interaction. The best energy res- 
olution and thus lowest energy threshold is obtained with 
the 4 fj,s shaped signal. A third branch of the energy out- 
put is shaped with 3 fj,s and amplified to record events 
up to 8 MeV in order to identify radioactive impurities 
contributing to the background. The spectra are mea- 
sured with 13bit ADCs, which also release the trigger 
for an event by a peak detect signal. Further triggers 
are vetoed until the complete event information, includ- 
ing the pulse shape, has been recorded. To record the 
pulse shape (for details see ]l3|] ) the timing output of 
the preamplifier is divided into four branches, each signal 
being integrated and differentiated in filtering amplifiers 
(TFAs) with different time constants. The TFAs are used 
since the charge current is integrated within the pream- 
plifier. The signals are amplified to record low as well as 
high-energetic pulses. 



DATA ANALYSIS 



An energy threshold of 9 keV has been reached. This 
rather high value is due to the large detector size and 
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a 50 cm distance between FET and detector. Both ef- 
fects lead to a higher system capacitance and thus to an 
enhancement of the baseline noise. 

We calibrate the detector with a standard 152 Eu- 228 Th 
source. The energy resolution at 727 keV is (2.37 ± 0.01) 
keV. In order to determine the energy resolution at 
keV the dependence of the full width at half maximum 
(FWHM) on the energy is approx imated with an em- 
pirical function y — y/a + bx + cx 2 (y = resolution, x 
= energy) in a \ 2 nt. The best fit (x 2 /DOF = 0.09) is 
obtained for the parameters a = 3.8, b = 2.2xl0 -3 , c = 
5xl0 -7 . The zero energy resolution is (2 ± 0.01) keV. 

To determine the energy threshold, a reliable energy 
calibration at low energies is required. The lowest ener- 
getic line observed in the detector was at 74.97 keV (K Q 
line of 208 Pb). This is due to the rather thick copper 
shield of the crystal (2 mm) which absorbs all low energy 
7 lines. Thus an extrapolation of the energy calibration 
to low energies is needed, which induces an error of (1- 
2) keV. Another possibility is to use a precision pulser 
in order to determine the channel of the energy spec- 
trum which corresponds to zero voltage pulser output. 
Since the slope of the energy calibration is independent 
of the intercept and can be determined with a calibra- 
tion source, this method yields an accurate value for the 
intercept of the calibration. The same method to deter- 
mine the offset of the calibration is also used by Q and 
Q. The pulser calibration reduces the extrapolated 9 
keV threshold systematically by (1-2) keV. In order to 
give conservative results, we use the 9 keV threshold for 
data analysis. 

In Fig. the energy deposition as a function of time 
is plotted. Accumulation of events (bursts) with energy 
depositions up to 30 keV can be seen. They are irregu- 
larly distributed in time with counting rate excesses up 
to a factor of five. These events can be generated by 
microphonics (small changes in the detector capacitance 
due to changes in temperature or mechanical vibrations) 
or by a raised electronic noise on the baseline. Although 
rare, they lead to an enhancement of the count rate in the 
low energy region. A possibility to deal with microphon- 
ics would be to exclude the few days with an enhanced 
count rate from the sum spectrum like in §. This would 
lead, however, to unwanted measuring time losses. An- 
other time filtering method is applied in the following. 

The complete measuring time is divided into 30- minute 
intervals and the probability for N events to occur in one 
time interval is computed for energy depositions between 
9-850 keV. In the histogram in Fig. the physical events 
are Poisson distributed. The mean value of the distribu- 
tion is (2.65 ± 0.06) events/30 min and cr=(1.67± 0.05). 
The cut is set at the value N = 2.65 + 3cr«8. With 
this cut less than 0.01% of the undisturbed 30 minutes 
intervals are rejected. The initial exposure of the mea- 
surement was 0.7 kgyr, after the time cut the exposure is 
0.69 kgyr. In this way, more than 98% of the initial data 
are used. A similar method to reduce the microphonic 



noise in the low energy region was also applied by || . 

Another way to reject microphonic events would be 
to analyze the recorded pulse shapes of each event. In 
a former experiment we showed [y_3| that pulse shapes 
of nuclear recoil and 7 interactions are indistinguish- 
able within the timing resolution of Ge detectors. Thus 
a reduction of 7-ray background based on pulse shape 
discrimination (PSD) is not possible. Consequently 7 
sources can be employed to calibrate a PSD method 
against microphonics, which was shown to reveal a differ- 
ent pattern in the pulse shape. Since such a pulse shape 
analysing method is still under development, we use only 
the Poisson-time-cut method in this paper. 

Figure shows the sum spectrum after the time cut. 
The background counting rate in the energy region be- 
tween 9 keV and 30 keV is 0.081 cts/(kgdkeV) [between 
15 keV and 40 keV: 0.042 cts/(kgdkeV)]. This is about a 
factor of two better than the background level reached by 
H with the same Ge detector. Table | gives the number 
of counts per 1 keV bin for the energy region between 
(9-50) keV. The dominating background contribution in 
the low-energy region from the U/Th natural decay chain 
can be identified in Fig. via the 352 keV and 609 keV 
lines (the continuous beta spectrum from 210 Bi originates 
from this chain). 



DARK MATTER LIMITS 



The evaluation for dark matter limits on the WIMP- 
nucleon cross section cr^T^ follows the conservative as- 
sumption that the whole experimental spectrum consists 
of WIMP events. Consequently, excess events from cal- 
culated WIMP spectra above the experimental spectrum 
in any energy region with a minimum width of the en- 
ergy resolution of the detector are forbidden (to a given 
confidence limit). 

The parameters used in the calculation of expected 
WIMP spectra are summarized in Table 0. We use 
formulas given in the extensive reviews |l5,|l6| for a 
truncated Maxwell velocity distribution in an isothermal 
WIMP-halo model (truncation at the escape velocity, 
compare also @). 

Since Ge is a spin zero nucleus, we give cross section 
limits for WIMP-nucleon scalar interactions only. For 
these, we used the Bessel form factor (see jl6| and refer- 
ences therein) for the parametrization of coherence loss, 
adopting a skin thickness of 1 fm. 

Another correction which has to be applied for a semi- 
conductor ionization detector is the ionization efficiency. 
There exist analytic expressions || for this efficiency, 
especially for germanium detectors and multiple experi- 
mental results measuring this quantity (see p3[ and ref- 
erences therein) . According to our measurements fl3| we 
give a simple relation between visible energies and re- 
coil energies: E v - ls = O.lAE^l^. This relation has been 
checked for consistency with the relation from B in the 
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relevant low energy region above our threshold. 

After calculating the WIMP spectrum for a given 
WIMP mass, the scalar cross section is the only free pa- 
rameter which is then used to fit the expected to the 
measured spectrum (see Fig. ) using a one-parameter 
maximum-likelihood fit algorithm. According to the un- 
derlying hypothesis (see above) we check during the fit 
for excess events above the experimental spectrum (for 
a one-sided 90% C.L.) using a sliding, variable energy 
window. The minimum width of this energy window is 
5 keV, corresponding to 2.5 times the FWHM of the de- 
tector (6tr width). The minimum of cross section val- 
ues obtained via these multiple fits of the expected to 
the measured spectrum gives the limit. Figure shows 
a comparison between the measured spectrum and the 
calculated WIMP spectrum for a 100 GeV WIMP mass. 
The solid curve represents the fitted WIMP spectrum us- 
ing a minimum width of 5 keV for the energy window. 
The minimum is found in the energy region between 15 
keV and 20 keV. The dashed line is the result of the 
fit if the energy window width equals the full spectrum 
width. It is easy to see that in this case the obtained 
limit would be much too conservative, leading to a loss 
of the information one gets from the measured spectrum. 

CONCLUSIONS 



actions are the most stringent limits obtained so far by 
using essentially raw data without background subtrac- 
tion. 

An improvement in sensitivity could be reached af- 
ter a longer measuring period. Higher statistics would 
allow the identification of the various radioactive back- 
ground sources and would open the possibility of a quan- 
titative and model-independent background description 
via a Monte Carlo simulation. Such a background model 
has already been established for three of the enriched Ge 
detectors in the Heidelberg-Moscow experiment and has 
been successfully applied in the evaluation of the 2i//3/3 
decay Q . A subtraction of background in the low en- 
ergy region in the form of a phenomenological straight 
line based on a quantitative background model for the 
full energy region (9 keV - 8 MeV) would lead to a further 
improvement in sensitivity. Background subtractions for 
dark matter evaluations of Ge experiments were already 
applied by |p3[ . 

Another way to reject radioactive background origi- 
nating from multiple scattered photons would be an ac- 
tive shielding in the immediate vicinity of the measuring 
crystal. This method is applied in our new dark mat- 
ter search experiment, HDMS (Heidelberg Dark Matter 
Search) pi] ]. HDMS is situated in the Gran Sasso Un- 
derground Laboratory and started operation this year. 



The new upper limit exclusion plot in the <r*^ r versus 
Mwimp plane is shown in Fig. . Since we do not use any 
background subtraction in our analysis, we consider our 
limit to be conservative. We are now sensitive to WIMP 
masses greater than 13 GeV and to cross sections as low 
as 1.12xl0~ 5 pb (for p = 0.3 GeV/cm 3 ). 

At the same time we start to enter the region 
(evidence-contour) allowed with 90% C.L. if the prelimi- 
nary analysis of 4549 kg days of data by the DAMA Nal 
Experiment [l§| ] are interpreted as an evidence for an an- 
nual modulation effect due to a spin independent coupled 
WIMP. 

Should the effect be confirmed with much higher statis- 
tics (20 000 kg days are now being analyzed by the 
DAMA Collaboration ]jl9|) it could become crucial to 
test the region using a different detector technique and a 
different target material. 

Also shown in the figure are recent limits from the 
CDMS Experiment from the DAMA Experiment 

p0| , as well as expectations for new dark matter exper- 
iments like CDMS Q, HDMS |ll and for our recently 
proposed experiment GENIUS |22j| . Not shown is the 
limit from the UKDM Experiment jllj which lies some- 
where between the two germanium limits. 

After a measuring period of 0.69 kgyr with one of the 
enriched germanium detectors of the Heidelberg-Moscow 
experiment, the background level decreased to 0.0419 
counts/ (kg d keV) in the low energy region. The WIMP 
nucleon cross section limits for spin-independent inter- 
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FIG. 1. Energy spectrum as a function of time for the enriched 76 Ge detector for the whole measurement period. Irregular 
bursts up to 30 keV can be seen. 
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FIG. 2. Distribution of the number of events with energies between 9 keV and 850 keV per 30-minute interval. The solid 
curve represents the Poisson fit to the data. 
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FIG. 3. Full recorded sum spectrum with lifetime 0.69 kgyr. Peaks from the U chain can be identified. A Cs peak 
slowly appears at 662 keV. All structures at low energies are fluctuations so far and therefore not identifiable. 



TABLE I. Number of counts per 1 keV energy bin after 250.836 kgd. 



Bin [keV] 


Counts [1/keV] 


Bin [keV] 


Counts [. 


9-10 


61 


30-31 


10 


10-11 


56 
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11 


11-12 


36 
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29 
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10 
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FIG. 4. Comparison of the measured low energy spectrum (shown from threshold to 50 keV) from the enriched 76 Ge Detector 
and calculated WIMP spectra for a 100 GeV WIMP, already fitted for the allowed cross section cr^^. The solid curve shows 
the result from the sliding-window fit (see text). The dashed curve would result for a full WIMP-spectrum fit to the data, 
yielding too conservative limits. 
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FIG. 5. Comparison of already achieved WIMP-nucleon scalar cross section limits (solid lines): the Heidelberg-Moscow 76 Ge 
M, the recent CDMS nat. Ge [go) and the new DAMA Nal result [juj, including their evidence contour [ jl8| , in pb for scalar 
interactions as function of the WIMP mass in GeV and of possible results from upcoming experiments (dashed lines for HDMS 
CDMS (at different locations) and GENIUS @). These experimental limits are also compared to expectations (scatter 
plot) for WIMP neutralinos calculated in the MSSM framework with non-universal scalar mass unification 
dots scaled for p = 0.5 GeV/cm 3 for comparison of published data). 
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TABLE II. List of parameters used for calculating WIMP spectra. 

Parameter Value 

WIMP velocity distribution 270 km/s 

Escape velocity 580 km/s 

Earth velocity 245 km/s 

WIMP local density 0.3 GeV/cm 3 

enriched 76 Ge mass 75.14 g/mol 
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